The link between the Hidden Broad Line Region and the accretion rate in 

Seyfert 2 galaxies 



Andrea Marinucci 1,2 , Stefano Bianchi 1 ' 3 ' 4 , Fabrizio Nicastro 2 ' 4 ' 5 , Giorgio Matt 1 , Andy D. 

Goulding 2 

Received / Accepted 

ABSTRACT 

In the past few years more and more pieces of evidence have been presented for 
a revision of the widely accepted Unified Model of Active Galactic Nuclei. A model 
based solely on orientation cannot explain all the observed phenomenology. In the 
following, we will present evidence that accretion rate is also a key parameter for the 
presence of Hidden Broad Line Regions in Seyfert 2 galaxies. Our sample consists of 
21 sources with polarized Hidden Broad Lines and 18 sources without Hidden Broad 
Lines. We use stellar velocity dispersions from several studies on the Call and Mg 
b triplets in Seyfert 2 galaxies, to estimate the mass of the central black holes via the 
Mbh-cfi, relation. The ratio between the bolometric luminosity, derived from the intrin- 
sic (i.e. unabsorbed) X-ray luminosity, and the Eddington luminosity is a measure of 
the rate at which matter accretes onto the central supermassive black hole. A separa- 
tion between Compton-thin HBLR and non-HBLR sources is clear, both in accretion 
rate (logL bo i/L Edd = -1.9) and in luminosity (logL bo i = 43.90). When, properly 
luminosity-corrected, Compton-thick sources are included, the separation between 
HBLR and non-HBLR is less sharp but no HBLR source falls below the Eddington 
ratio threshold. We speculate that non-HBLR Compton-thick sources with accretion 
rate higher than the threshold, do possess a BLR, but something, probably related to 
their heavy absorption, is preventing us from observing it even in polarized light. 

Our results for Compton-thin sources support theoretical expectations. In a model 
presented by Nicastro (2000), the presence of broad emission lines is intrinsically 
connected with disk instabilities occuring in proximity of a transition radius, which is 
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a function of the accretion rate, becoming smaller than the innermost stable orbit for 
very low accretion rates and therefore luminosities. 

Subject headings: Galaxies: active - Galaxies: Seyfert - Galaxies: accretion 



Introduction 



The widely accepted unification model f or active galactic nuclei (AGNs) invokes the same 
paradigm for Seyfert 2 and Seyfert 1 galaxies (|Antonuccilll993|) . The two different types of galax- 
ies are believed to be intrinsically the same but differ, observationally, due only to orientation: 
Seyfert 1 galaxies are observed at large accretion disk-line of sight angles (face-on) while Seyfert 
2 galaxies are seen edge-on through large column densities of obscuring material in the accretion 
disk plane, which prevent the direct view of the nuclear regions of these sources. This scenario 
came into existence after the discovery of polarized broad permitted emission lines in one of th e 
brightest Seyfert 2 galaxies, NGC 1068 (lAntonucci and Miller!! 1985k Miller and Antonuccill 1983k 
suggesting a geometry in which: (a) BLRs are confined in a relatively small region (~ light-days) 
surrounding the central source, (b) their direct view (in Seyfert 2s) is obscured by a flat distribution 
of distant material coplanar with the disk-plane, and (c) their line emission is Compton- scattered 
into the line of sight direction off a population of hot electrons extending at large radii above and 
below the accretion disk. 

In the past few years, exceptions to orientation-based unification models have been found, 
suggesting the possibilty that not all Seyfert 2 galaxies host a Seyfert 1 nucleus. Indeed, spec- 
tropolarimetric surveys find that only about the 50% of the bright est Seyfert 2 gala xies show hidden 
broad-line regions (HBLR) in their optical-polarimetric spectra dTranlboO lL Eooih . Several authors 
have suggested that the presence of BLRs in Seyfert 2s can be linked to the luminosi t y of the ac 
tive nucleus, and may disappear at low luminos i ties ( 



Gu and Huangll2002k Ma rtocchia and Matt 2002 



Tran 



.umsden and Alexander 2001 



2003 



Laor 



2003 



Tran 



Elit zurand Hoi 12009 ) or 



2001; 



low accretion rates (|Nicastrol 120001 : iNicastro et al.l 12003k ICzerny et al.l l2004f) . In both cases, the 
presence of the BLR is not an inescapable feature of all Seyfert galaxies, as postulated by the uni- 
fication model, but it is tightly linked to a physical parameter of the AGN, either the l uminosity or 
the accretion rate. However, recent works, casted doubts on these conclusions (e.g. iBian and Gu 
20071) . In these works large, but inhomogeneous, samples were used to identify a physical pa- 
rameter responsible for the existence or absence of HBLR in AGNs, but no clear indication for 
the existence of such a parameter was found. We think that this is largely due to the use of the 
[OIII] luminosities as a proxy of the nuclear activity of the AGN, and the difficulty in correcting 
this for extinction in obscured objects. It has been shown by several author, by comparing AGN 
luminosities derived from [OIII] and [OIV] emission line luminosities, that the observed [OIII] 
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luminosities often suffer significant attenuation in typical Seyfert 2 galaxies (e.g.lHaas et al 



Melendez et a 



Kraemer et al 



20081 : biamond-Stanic et al.ll2009l : boulding & Alexanderll2009l : lBaum et al 



2005 



20101 : 



201 II) and that the efforts to use standard extinction co rrection for [OIII] are not al 



ways reliable (e.g. lGoulding & Alexanderj|2009l : LLaMassa et al.ll2010|) . 



[OIV] emission line luminosities might therefore be more reliable than [OIII] luminosities, as 
tracers of the intrinsic nuclear emission. However, they are still an indirect proxy of the nuclear 
continuum emission. In this paper, instead, we use the observed X-ray (2-10 keV) continuum 
emission as a direct probe of the AGN activity. We selected a sample of type 2 AGN with good 
quality spectropolarimetric and X-ray observations, for which we can give a good estimate of the 
mass of the central s upermas si ve black hole (BH).[jBy doing so, we find evidence suggesting that 
accretion rate is the main parameter that sets the existence of HBLRs in Seyfert 2 galaxies. 



2. The Sample 



O ur starting sample is mostly based on the spectropolarime tric surveys performed by iTran 



d 1 995L 1200 lL l2003h on the Seyfert 2 galaxies included in the CfA (IHuchra and Burgll 19921) and 12 
fim ( Rush et all 1993 ) samples. Additionally, we included objects from other, high-quality, spec- 
tropo larimetric studies ( Tran et a0l992 : Young et aDll996 : Moran et al. I2OO0I : Lumsden and Alexander 
200 II) . From this sample, consisting of 90 candidates, we selected only sources with available mass 
(Section 2.1) and bolometric luminosity (Section 2.2) estimates, and so for which their accretion 
rate can be evaluated. For each source, we evaluate their accretion rate in units of Eddington, as 
the ratio L hol /L Edd (hereafter A Edd ), where 



L Edd = 1.2 x 10 



38 



M 



BH 



Mr 



erg s 



(1) 



2.1. The BH mass sub-selection 



In this work, B H masses are homoge nuously derived for the entire sample, by using uniquely 
the Msn-cr* relation (|Tremaine et al.ll2002l) : 



M BH = 1.35xl0 8 ( 



CT* \4.02 



200 km r 



(2) 



'Our work is based on a previous study discussed in lNicastro et al.l (12.003b . but here we use a wider sample and 
more robust estimates of the BH masses. 



Table 1 : The sample used in this paper. X-ray and bolometric luminosities are already corrected for a factor 70 for all Compton Thick sources. Columns: 
(a) Name of the object; (b) Instrument (AS: Chandra ACIS-S; PN: XMM-Newton EPIC pn; XIS: Suzaku XIS0+3) and OBSID; (c) /if Compton-thin, X 



if Compton-thick; (d) Stellar velocity dispersions; (e) Mass of the central BH in solar mass units; (f) Absorbing column density in 10 cm units (g) 
Absorption corrected 2-10 keV luminosity; (h) Bolometric luminosity; (i) Accretion rate in Eddington units; (1) Extinction-cor rected O full (A 5007 ) 
lumino sity; (m) O fiv] (/125.89 um ) lu min osity (n) Reference for stell ar velocity d i spersio ns, HBLR (in order); R efer ences: (1) 



(2005), (2) 



Nelson and Whittle 
19951) . r81 lTranl(l2003l) . (9 faanl 



19951) (3) lTerlevichetlii1(ll990h . (4) 
d 1995b . (TOl llvloran etalTlhoOOl) . (11) 



Oliva et al 



Garcia-Rissmann et al 



Jl999h (5) iBarth et all (booi. (6) Tshaw et all dl993h. (7) iMcElrov 



Young et al.1 (119961) . (12) lLumsden and Alexander! (1200 II) . (13) Other surveys 



Hidden Broad Line Region Seyfert 2s 

Object name Obs. ID C-Thin cr, (kms -1 ) log(M w ,/M ) N H log(L 2 -io kcv) \og(L Bo i) log(A edd ) log(L [0 ///]) log(L [0/V /j) References 

(a) (b) (c) (d) (e) (f) (g) (h) (i) (1) (m) (n) 



CIRCINUS 


PN- 0111240101 


X 


75 ±5 


6.42 


430+*° 

— 


42.62 


43.76 


-0.75 


40.92 


40.58 


1, 13 


IC 3639 


XIS-702011010 


X 


99 ±5 


6.90 


>150 


42.64 


43.80 


-1.20 


41.89 


40.66 


1,8 


IC 5063 


XIS-704010010 


/ 


160 ± 25 


7.74 


25+3 


42.83 


44.03 


-1.81 


41.56 


41.40 


2.8 


IRAS01475-0740 


PN-0200431101 


X 


108 ± 17 


7.05 


> 200 


43.52 


44.90 


-0.25 


41.76 


40.74 


1.8 


MCG -2-8-39 


PN-0301 150201 


y 


126 ± 11 


7.32 


44+g 


42.57 


43.70 


-1.72 


41.16 


41.07 


1.8 


Mrk3 


XIS-100040010 


X 


249 ±4 


8.51 


140+?, 


44.44 


46.10 


-0.51 


43.27 


41.97 


1, 10 


Mrk 348 


PN-0067540201 


✓ 


117 ± 18 


7.19 


15!? 


43.41 


44.77 


-0.52 


41.96 


40.95 


2,8 


Mrk 1210 


PN-0002940701 




82 ± 16 


6.57 


22% 


43.02 


44.24 


-0.43 


42.37 




1.9 


NGC513 


PN-0301 150401 


✓ 


150 ± 25 


7.63 




42.72 


43.93 


-1.80 


41.14 


40.74 


2,8 


NGC 591 


PN-0200431001 


X 


107 ± 18 


7.04 


> 160 


43.02 


44.27 


-0.87 


41.97 




2,10 


NGC 788 


XIS-703032010 


/ 


140 ± 20 


7.51 


85 -6 


43.15 


44.43 


-1.18 


40.73 


40.97 


2,10 


NGC 1068 


PN-01 11200201 


X 


147 ±3 


7.59 


> 1000 


43.02 


44.27 


-1.42 


42.38 


41.81 


1,8 


NGC 2273 


XIS-702003010 


X 


136 ± 22 


7.46 


120+>>° 


42.73 


43.90 


-1.65 


41.13 


40.07 


2, 10 


NGC 3081 


XIS-70301301 


/ 


113 ±4 


7.13 


° y -l5 


42.50 


43.61 


-1.62 


41.43 


41.09 


1,10 


NGC 4388 


XIS-800017010 


/ 


119 


7.22 


31+| 


42.90 


44.12 


-1.20 


41.85 


41.58 


3,8 


NGC 4507 


XIS-702048010 


✓ 


152 ±4 


7.65 


87 -8 


43.11 


44.39 


-1.36 


42.19 


41.02 


1,10 


NGC 5252 


PN-0152940101 


✓ 


190 ± 27 


8.04 


z - z -o.i 


43.04 


44.30 


-1.84 


42.05 




2,11 


NGC 5506 


PN-0554170101 


✓ 


180 


7.95 


3 0+ 01 


43.05 


44.30 


-1.75 


41.45 


41.28 


4,8 


NGC 7212 


PN-0200430201 


X 


140 ±9 


7.51 


> 150 


43.77 


45.22 


-0.38 


42.73 




1,13 


NGC 7674 


PN-0200660101 


X 


144 ± 32 


7.56 


> 100 


43.96 


45.47 


-0.18 


42.57 


41.97 


2,8 


NGC 7682 


PN-0301 150501 


X 


123 ± 17 


7.28 


> 100 


43.02 


44.27 


-1.11 


41.76 


41.01 


2.8 



Non Hidden Broad Line Region Seyfert 2s 



Object name 


Obs. ID 


C-Thin 


<t„ (km s ) 


log(M w ,/M Q ) 


N H 


log(L 2 -l0kev) 


log(L Bo ,) 


log(/W) 


log(i[o/;/]) 


l°g(£[0/V]) 


References 


(a) 


0) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


(i) 


(1) 


(m) 


(n) 


M51 


PN-0303420101 


X 


82 ± 11 


6.57 


>400 


41.54 


42.51 


-2.17 


40.03 


39.01 


2,8 


Mrk 573 


AS-7745 


X 


148 ± 3 


7.60 


> 100 


43.10 


44.37 


-1.33 


42.39 


41.71 


1.8 


Mrk 1066 


PN-0201770201 


X 


1 19 ± 19 


7.22 


> 100 


42.92 


44.15 


-1.17 


42.27 




2,10 


NGC 1320 


PN-0405240201 


X 


124 ± 14 


7.29 


> 100 


42.69 


43.86 


-1.54 


41.08 


40.67 


2,8 


NGC 1358 


PN-0301650201 


X 


185 ±20 


7.99 


, 30+850 


43.05 


44.31 


-1.79 


41.36 


40.56 


2,10 


NGC 1386 


XIS-702002010 


X 


133 ±3 


7.42 


> 100 


41.62 


42.59 


-2.92 


41.09 


40.25 


1,8 


NGC 1667 


XIS-701006010 


X 


173 


7.88 


> 100 


42.56 


43.70 


-2.28 


42.03 


40.51 


3,8 


NGC 3079 


XIS-803039020 


X 


150 ± 10 


7.63 


> 100 


42.02 


43.05 


-2.68 


40.48 


39.53 


6.8 


NCG 3281 


XIS-703033010 


✓ 


176 ±3 


7.91 


79 +io 


42.65 


43.81 


-2.20 


41.30 


41.58 


1,10 


NGC 3982 


PN-0204651201 


X 


81 ± 13 


6.55 


> 100 


41.15 


42.06 


-2.59 


40.33 


39.10 


5,8 


NGC 5135 


XIS-702005010 


X 


124 ±6 


7.29 


> 100 


43.10 


44.37 


-1.02 


42.21 


41.48 


1,8 


NGC 5283 


AS-4846 


/ 


148 ± 14 


7.60 


12+jj 


41.54 


42.49 


-3.21 


40.88 




2,8 


NGC 5347 


XIS-703011010 


X 


103 ± 14 


6.97 


> 100 


42.39 


43.49 


-1.58 


41.22 


40.01 


3,8 


NGC 5728 


XIS-701079010 


X 


209 


8.21 


> 150 


43.29 


44.61 


-1.69 


42.83 


41.36 


7,11 


NGC 7582 


PN-01 12310201 


/ 


1 13 ± 3 


7.13 


57+ 10 


41.90 


42.91 


-2.32 


41.63 


41.11 


1,8 


NGC 7130 


XIS-703012010 


X 


147 ±5 


7.59 


> 100 


43.10 


44.37 


-1.32 


42.55 


40.94 


1,12 


NGC 7172 


PN-0202860101 


/ 


160 ±9 


7.74 


8 4+0.1 
°- -0.1 


42.73 


43.90 


-1.93 


40.84 


40.82 


1,8 


UGC 6100 


PN-0301 15 1101 


X 


156 ± 25 


7.69 


> 100 


42.84 


44.05 


-1.75 


42.18 




2,8 
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Stellar velocity dispersions (<x*) were mainly taken from LNelson and Whittle! (1 1995|) and lGarcia-Rissmann et 



(|2005|) . These estimates are all based on direct meaurements of Can (8498 A, 8542 A and 8662 

o □ o o 

A) and Mgl 3p — > 4s (5167 A, 5172 A and 5183 A, hereinafter Mg b) triplet absorption, imprinted 
by the interstellar medium of Seyfert 2 galaxies. For sources with more than one measurement, we 
used that with smaller error bars. Uncertainties on the BH mass estimates based on the M B h-o~* 
relation, come from the statistical errors on the cr* measureme nts (listed in Table 1) as well as from 
the spread in the Mbh-ct* relation itself, estimated in 0.44 dex (Giilte kin et al. 2009). This spread is 



generally much larger than the statistical error on ov We therefore assumed an uncertainty of 0.44 
dex for all our BH mass estimates. Stellar velocity dispersions and their associated uncertainties, 
as well as BH masses, are listed in Table \T\ (columns (d) and (e), respectively). 

The cr* selection reduced the original spectropolarimetric sample to 46 sources with direct cr t 
measurements. 



2.2. The X-Ray Sub- Selection and the Final Sample 

To estimate the bolometric luminosity of the sources of our sample, we use the 2-10 keV 
luminosity, which is a direct tracer of the primary emission, and apply a bolometric correction. We 
searched the XMM-Newton, Chandra, Suzaku (the front-illuminated CCDs have the largest area 
at high energies), and Swift archives for observations of the sources of our sample of 46 objects 
with accurate BH-mass estimate. When multiple observations of a single target were available, we 
selected that with the highest S/N ratio in the 2-10 keV band. For all sources of our final sample, 
a minimum of 150 counts in the 2-10 keV band was required, to derive reliable estimates of the 
column density and, therefore, intrinsic X-ray luminosity. These criteria reduced the final sample 
to a total of 39 sources: 21 with polarized hidden Broad Emission Lines (HBLRs) and 18 sources 
without (non-HBLRs). These are listed in Tabled] together with the selected X-ray observations. 



3. Data Analysis 

3.1. X-ray Data Reduction and Analysis 

Suzaku X-ray Imaging Spectrometer (XIS) data were processed with the latest calibration files 
available at the time of the analysis (2011-02-10 release), by using FTOOLS 16 and SUZAKU 
software Version 2.3, and adopting standard filtering procedures. Response matrices and ancil- 
lary response files were generated using XISRMFGEN and XISSIMARFGEN. The 0.5-10 keV 
spectra extracted from the front-illuminated XIS0 and XIS3 have been co-added via the ftool AD- 
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DASCASPEC. 



We used only XMM-Newton observations performed with the EPIC-Pn camera (|Striider et al 



200 1|) operated in large window and medium filter modes. Source data 'cleaning' (exclusion of 
flaring particle back ground intervals) and spectra extraction, were performed with SAS 10.0.0 
([Gabriel et al.ll2004l) via an iterative proces s which leads to a maxi mization of the Signal-to-Noise 
Ratio (SNR), similarly to that described in IPiconcelli et al.l (120041) . For each source, background 
spectra were extracted from source-free circular regions of the source field. 

Fi nally, Chandra data w ere reduced with the Chandra Interactive Analysis of Observations 
(CIAO: lFruscione et alJ2006h 4.3 and the Chandra Calibration Data Base (CALDB) 4.4. 1 software, 
by adopting standard procedures. 

All the spectra with a high S/N ratio were binned in order to over- sample the instrumental 
resolution by at least a factor of 3 and to have no less than 30 counts in each background-subtracted 
spectral channel. This allows the applicability of the x 1 statistics. The adopted cosmological 
parameters are H = 70 km s~ ! Mpc" 1 , Q. A = 0.73 and Q m = 0.27. Errors are quoted at a 
confidence level of 90% for one interesting parameter (A^ 2 = 2. 7), if not othe rwise stated. The 
spectral analysis was performed with the package XSPEC 12.7.0 (|Arnaudlll996|) . 



For all sources of our sample, we fit their 2-10 keV spectra with a general baseline model, 
consisting of a power law continnuum atennuated by the line of sight column of Galactic absorbt- 
pion plus intrinsic absorption at the source redshift, plus three additional emission components: (a) 
a photoionized plasma emitter, to model the soft excess often detected in Seyfert 2s at E< 2 keV, 
(b) a cold Compton Reflector, scattering the primary nuclear photons off the inner walls of cold 
circum-nuclear material and along the line of sight, and (c) positive gaussian profiles, to model 
fluorescence emission lines of high-Z elements such as Fe at 6.4 keV, as required by the data. The 
model can be parameterized as: 



F(E) 



e'^ E)N hPh c + e 



-cr(E)N H 



BE~ T + R(T) + ^ Gi(E)] 



(3) 



where cr(E) is the photoelectric cross-sectio n (abundances as inlAnders and Grevessdl 198 9"). N% 
is the line-of-sight Ga lactic column densit y (IDickey and Lockmanll 19901) : Ph c is the photoionized 
plasma emission (see LBianchi et al.l |2010L for details on the adopted CLOUDY model); Nh is the 
neutral absorbing column density at the redshift of the source; B is the normalization of the pri- 
mary po werlaw with slope Y; R(T) is the Compton-reflection component (modelled in XSPEC with 
pexrav (|Magdziarz and Zdziarskilll995l) ): and Gi(E) are the required Gaussian profiles. 



Our adopted models do not aim at obtaining the 'best fit' for each source, but only a reliable 
estimate of the observed 2-10 keV luminosity (Tabled] column (g)). A more detailed description 
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of the results of the fits, in particular for the photoionization CLOUDY component, is deferred to a 
future work. 



3.2. Nuclear 2-10 keV Luminosity for Comp ton-Thick Sources 



Compton-Thick sources are defined as sources for which the intrinsic column of absorbing 
gas, Nh, is greater than 10 24 cm -2 , so preventing us from directly observing the primary pow- 
erlaw below 10 keV. The 2-10 keV spectra of Compton-thick sources are therefore generally ade- 
quately modeled by a pure-reflection component (i.e. reflection-dominated Seyfert 2s), which can 
only provide a measurement of the reflected, and not the nuclear, luminosity. 
Severa l tracers can be used to estimate t he bolometric lu minosity of Compton-thick objects, e.g. 
[Om] (ILamastra et al.ll2009l) and [OlV] (IRigby et al.ll2009l) . For homogeneity with the Compton- 
thin sources of our sample, here we use the observed X-ray luminosity, by applying the following 
procedure. 

We compare the extincti on-corrected [Ol lll l uminosities of the s ources of our sample (Ta- 
ble[D column (1)), taken from lWu et al.1 (|201l|) and lBian and Gul (|2007|) . with their observed 2-10 
keV reflection-dominate d luminosities. Bot h works correct the observed F[ ni] for extinction, by 
applying the relation in lBassani et al.l (|1999|) : 



pcor 
r [OIII] 



.[ 



(Ha/H0) obs ,2.94 



[0 '"H (HalHftk 



(4) 



where an intrinsic Balmer decrement (Ha/Hf3) =3 is adopted. The results are shown in Fig.[H top 
panel. A ll the Compton-Thick sources have a L x /L[oni] rat i° smaller than the best fitting relation 
found by ILamastra et al.l (|2Q09r) for a large sample of Compton-thin Seyfert 2s (solid line in Fig.Q] 
compare with the dotted line in the same figure, which best-fits the sub-sample of Compton-thin 
sources of our sample). For the Compton-thick sources of our sample, we find a mean value 
of < log^y/L^TTTj) >= -0.76 ± 0.09. For the Compton-thin sources of their sample, instead, 
Lamastra et al.1 (I2009|) found a mean < log(Lx/L[oiii]) >= 1-09. The ratio between these two 
means can be used as an estimate of the correction factor needed to infer the nuclear 2-10 keV 
luminosity of Compton-thick sources, from their observed 2-10 keV luminosity: 1 o 109 ~ ( ~°- 76) ~ 70. 
The bottom panel of Fig. Q] shows the logL x - /ogL [0 ni] plane for all the sources of our sample, 
after applying a correction factor of 70 to the observ ed 2-10 keV luminosit y of all Compton-thick 
sources. Now all the sources of our sample lie on the ILamastra et al.1 (I2009|) relation for Compton- 
thin sources. 



We decided to adopt this correction for all Compton-thick sources, also for those where E>10 
keV data were available (Suzaku PIN) and so a direct estimate of the nuclear continuum above 10 
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keV could in principle be attempted. This is because Compton scattering at high column densities 
may significantly suppress the observed intrinsic luminosity even at E> 10 keV, and the model- 
lization of these effects are highly dependent on the (unknown) geometry of the absorber (e.g. 



Matt et al.l 1 19990 . In the cases where objects presented observations in both Compton-thin and 



Compton-thick states in the past few years we chose the former, to infer with a higher precision 
the unabsorbed luminosity of the source. In one case (IRAS01475 -0740 ), the 2-10 keV source 
spectrum is ambiguous: following the detailed analysis presented in lBrightman and Nandral (I2008|) 
we classify this object as Compton-thick. 



3.3. Bolometric Luminosities and Eddington Ratios 



For each source of our sample, we used their intrinsic 2-10 keV luminosity (observed, for 
Compton-Thin sources, or inferred, for Compton-Thick sources, see previous section) to derive 
their bolometric luminosity, and thus Eddington ratio (Table [H colum ns (h) and (i)), by adopting 
the luminosity-dependent relation presented in iMarconi et al.1 (|2004|) . We tested the importance 



of this assumption in our calculatio ns, by trying different bolometric corrections (e.g. Elvis et al. 



19941 : IVasudevan and Fabianll2009|) . and comparing the results. The most significant differences, 
between these methods, were found at the extremes of the luminosity (accretion) range spanned 
by the sources of our sample. Therefore, the exact choice of the bolometric correction does not 
significantly affect our results. 

The major contribution to the bolometric luminosity error determination comes from the un- 
certainty on the bolometric correction itself, which is generally signific antly larger than the u n- 
certainty on the X-ray luminosity. The bolometric correction proposed by LMarconi et al.1 (|2004l) is 
based on the correlation between the UV luminosity at 2500 A a nd the 2-10 keV X-ray luminosity, 
whose spread is estimated as 0.37 dex lYoung et al.l (see eq. 3 of l2010h . Therefore, we assume this 
uncertainty on our derived bolometric luminosities. 

For the Eddington ratios, we propagated the uncertainties on the BH mass and bolometric 
luminosity. This leads to an uncertainty of 0.5 dex in the Eddington ratios. 



4. Results 

We first considered only the Compton-thin sources of our sample. In the top panel of Fig. [2l 
we plot the bolometric luminosity against the Eddington ratio for all the Compton-thin sources 
of our sample. HBLR and non-HBLR sources are clearly separated both in luminosity and in 
accretion rate. To address the statistical significance of this separation (both in luminosity and 
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Table 2: Results from the two-sample K-S test. For Compton Thick sources a correction factor of 
70 is used for the calculation of the 2-10 keV luminosity and hence for the bolometric luminosity. 
(1) Description of the sample; (2) number of objects in the sample; (3) percentage probability that 
the two samples derive from the same parent population in L Boi ; (4) percentage probability that the 
two samples derive from the same parent population in L 2 _io kev, L [0im , or L [om ; (5) percentage 
probability that the two samples derive from the same parent population in A edd ; (6),(7),(8) maximal 
separation value between the cumulative distributions of the two samples. 



Sample 


N 


P(L Bo i) 


P(L2-I0kev) 




i max 


j max 


A edd 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Non-HBLR Compton Thin 
HBLR Compton Thin 


4 
11 


0.5% 


4.4% 


0.1% 


43.90 


42.73 


-1.9 


Non-HBLR Compton Thick and Thin 


18 


1.8% 


6.9% 


1.0% 


43.90 


42.56 


-1.8 


HBLR Compton Thick and Thin 


21 


[OIII] 












j max 




Non-HBLR Compton Thin 
HBLR Compton Thin 


4 
11 


38.5% 


38.5% 


8.7% 


43.51 


41.36 


-2.5 


Non-HBLR Compton Thick and Thin 
HBLR Compton Thick and Thin 


18 
21 


23.2% 


23.2% 


25.7% 


43.54 


41.39 


-2.1 


[OIV] 






P(Wi) 






jmax 




Non-HBLR Compton Thin 


3 














HBLR Compton Thin 


9 














Non-HBLR Compton Thick and Thin 


15 


10.4% 


10.4% 


1.1% 


43.82 


40.57 


-1.4 


HBLR Compton Thick and Thin 


17 
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accretion rate), we performed a two-sample Kolmogorov-Smirnov test on our data. These tests 
give probabilities of 0.5% (luminosity) and 0.1% (Eddington ratio) that the two classes are taken 
from the same parent population. These correspond to 2-sided Gaussian-equivalent significances 
of 2.8<t and 3.3cr, respectively. The threshold Eddington ratio and luminosity that minimize the 
probability of chance separation, are log A e dd = -1.9 and logL bo i = 43.90, respectively. 

We then repeated our analysis, by including also Compton-thick sources (for which inferred 
intrinsic 2-10 keV luminosity estimates are uncertain: see Sect. 12.2b ■ The inclusion of Compton- 
thick sources makes the HBLR versus non-HBLR separation less sharp (Fig. [21 bottom panel), 
and less significant, with probabilities of chance separation of 1.0% and 1.8%, for Eddington ratio 
and bolometric luminosity, respectively (Table©, corresponding to statistical significances of only 
2.4cr and 2.6a. However, we note that no HBLR source falls below the Eddington ratio threshold. 



4.1. Comparison with other Works and Methodologies 

Other authors investigated larger (but, in most cases inhomogeneous) samples to search for 
physical parameters responsible for the existence or absence of HBLRs in AGNs , and found often 



no clear indication for the existence of such a parameter (e.g. iBian and Gull2007|) . However, we 
think that this is due to the non-homogenous nature of the assembled samples, or to the ill-suited 
methodologies used to derive bolometric luminosities and/or BH masses, or to a combination of 
both issues. 

BH masses can be derived through a variety of different means (e.g. FWHM[OIII], cr*). Using 
all of these possible means, obviously allows for large samples to be assembled. However, building 
samples with black- hole mass determinations derived with different methods, introduces signifi- 
cant sources of uncertainties, which may wash out possible correlations among important physical 
parameters, like the source accretion rate or luminosity. For example, stellar vel ocity dispersion 



can in principle be inferred through the relation cr* = (FWHM [0 ///] /2.35) x 1 / 1 .34 ([Greene and Ho 



2005), by making use of [Om] FWHM measurements, which are available for a large number 
of sources in the literaure. We tested the goodness of this estimator against the sources of our 
sample, for which direct stellar velocity dispersions are available. For all sour ces of our sam 



pie, we searched the literature for [Om] FWHM measurements, and applied the IGreene and Ho 



(|2005|) relation to derive the stellar velocity dispersions. We then compared these estimates with 
the measurements of cr*. Fig. [3] shows the result of this comparison: clearly the two values differ 
significantly and do not appear to be linearly correlated with one another. We, therefore, use only 
BH estimates based on direct measures of the cr*, in order to avoid further uncertainties due to the 
indirect derivation of this parameter. 
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Another important source of potential error, is the proxy used to derive the nuclear bolometric 
luminosity of an obscured AGN. In our work, we use 2-10 keV luminosities, which are a direct 
probe of the intrinisic nuclear activity. However, such luminosities are available only for a lim- 
ited number of Seyfert 2s. On the contrary [Olll] luminosities only echoes the intrinsic nuclear 
activity, but are available for a much larger number of obscured AGNs. To check the goodness 
of the [Olll] luminosity estimator, we compared our results with those obtained, for the sources 
of our sample, by using the the [Olll] luminosities and the b olometric luminosities derived from 
them, through the bolometric correction C[o///] described in lLamastra et al.1 (|2009f) . The results 
are shown in Table [2] and Fig @] (top panel): HBLR and non-HBLR sources are now mixed in 
both bolometric luminosities and Eddington ratios, with no significant separation between the two 
classes. The presence or absence of a clear separation, either in acrretion rates or bolometric lumi- 
nosities, between Compton-thin HBLRs and non-HBLRs, when these quantities are inferred from 
2-10 keV (probabilities of chance separation of 0.1% and 0.5%, for \ Edd and L Boh respectively) or 
rOllll (probabilities of chance separation of 8.7% and 38.5%) observed luminosities, confirms (e.g. 



Lamastra et al.ll2009l : iTrouille and BargejfeOlOL and references therein) that the [Olll] luminosity 



is not a direct proxy for the bolometric luminosity of an AGN. 

A large number of sources in our sample (32 out of 39) also have [Oiv] measurements (|Pereira-Santaella et al 



2010; Weaver et al 



2010L For a few s ource s we inferred [Oiv] lumi nosities following the m ethod 



presented in iGoulding & Alexander ( 2009 ). Using the relation in Goulding et al. ( 2010 ). both 
bolometric luminosities and accretion rates can be inferred from the [Oiv] luminosities. The re- 
sults are shown in the bottom panel of Fig. |4]and Table |2] The two classes show a more significant 
separation in accretion rate than in bolometric luminosity (probability of chance separation of 1 . 1 % 
and 10.4%, respectiveljH). Both HBLR and non-HBLR Compton thick sources are included in the 
sample and the separation between the two different populations in accretion rate resembles the one 
inferred from our X-ray analysis, sug gesting that the efforts to use standard extinction c orrection 
for [Oin] are not always reliable (e.g. IGoulding & Alexandern2009l: LLaMassa et al.ll2010|) . 



5. Discussion 



Since the seminal results by lMiller and Antonuccil(|l983|) on the archetypical Seyfert 2 galaxy, 
NGC1068, many efforts have been spent to understand whether the Unification Model is valid for 
all Seyfert galaxies, or the presence/absence of the BLR may be also linked to intrinsic pro pert ies of 
diffe rent classes of sources. Intrinsic differences were strongly advocated by lTranl(|200l|) and lTran 
(|2003l) . who found that HBLR Seyfert 2s are normally associated with typical obscured Seyfert 1 



2 A K-S test for the subsample of only C ompton-thin so urces cannot be performed, since the minimum number of 
data values for a statistical significance is 4 dstephensll97d pag. 120-121) 
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nuclei, while non-HBLR Seyfert 2 galaxies host, on average, significantly weaker muclei, likely 
incapable of generating classical BLRs. 

This evidence supports theoretical models that link the prese nce or absenc e of BLRs in 
Seyferts, to intrinsic nuclear properties. In a model proposed by iNicastrol (|2000l) . the presence 
of broad emission lines is intrinsically connected with disk instabilities occuring in proximity of 
a transition radius at which the accretion disk changes from gas-pressure dominated to radiation- 
pressure dominated. The transition radius is function of the accretion rate, and becomes smaller 
than the innermost stable orbit for accretion rates (and therefore luminosities) lower than a thresh- 
old that depends weakly on the BH mass. Weak AGN should, therefore, lack the BLR. 

In this work we try to test this model in the least biased observational way, by collecting the 
'cleanest' possible sample of Seyfert 2 galaxies with the best spectropolarimetric data available at 
the moment, and robust estimates of their BH masses performed homogeneously on the basis of the 
observed stellar velocity dispersion. When only Compton-thin sources are considered, the estimate 
on their accretion rate is done directly on the observed nuclear X-ray emission, and the modest ab- 
sorption is likely not able to affect severely the detection of the hidden BLR, if present. We find that 
the separation between HBLR and non-HBLR Compton-thin sources is highly significant both in 
accretion rate and luminosity. In particular, no HBLR is found below the threshold Eddington rate 
log A e dd = -1.9, and no non-HBLR above the same limit. Even when Compton-thick sources are 
included, no HBLR is still found at accretion rates lower than the above threshold. This threshold 



accretion rate is in good agreement with the value presented in a recent work (ITrump et al.ll201 ll) 



This result supports the theoretical expectations of the model proposed by INicastrol (120 00) 



albeit with a slightly higher value of the threshold accretion rate (log A e dd ~ -2.5 in 



Nicastro 



( 2000 )). We stress here that the model proposed by Nicastro (|200o|) applies o nly to radiatively ef- 
ficien t AGNs, accreting through a Shakura-Sunyaev disk (SS-disk hereinafter. IShakura & Sunyaev 
1973b . Bro ad Emission lines ar e known to exist in some objects (mostly LINERs o r transitio n 
objects, e.g. Elitzur and Ho 2009|) accreting at rates as low as log A ed d ~ -6 (e.g. M 81, Ho 2008 ). 
These objects are highly radiatevely inefficient, and therefore most likely do not host a classic 
SS-disk. 

If BLRs do not exist in weakly accreting AGN, one would expect the existence of unabsorbed, 
genuine Seyfert 2 galaxies. Such objects do indeed exist, and the best examples (where the lack of 
optical broad lines and of X-ray obscuration are unambiguosly assessed in simultaneous high S/N 
observations) ha ve Eddington rates lower than the threshold esti mated in this paper : NGC 3147 
(\ogA edd - -4: lBianchi et alJl2008ah . Q2131 427 (log/L^ -2.6:IPanessa et al. 
(logA e dd - -2: Bianchi et al., in preparation; Brightman and Nandral 2008 ). 



2009), NGC 3660 



It is interesting to note that among all the 156 X-ray unobscured (N H < 2 x 10 cm 2 ) AGN ob 
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served with XMM-Newton (CAIXA: Bianchi et al.ll2009 ). only 6 have an Eddington rate lower than 
the threshol d found in this worljjj. Four of these 6 sources were previously classified as una bsorbed 



Sev fert 2s (IPanessa and Bassani 



20021 : IPanessa et al J 12009 : Brightman and Nandral 12008) but see 



also lShi et al.1 (e.g. l2010l and references therein), who recently claimed the presence of very broad 
and weak broad emission lines in two of these objects). The remaining two objects are PG 101 1- 
040 and NGC 7213, and both show the presence of BLR in their optical spectra. Both objects are 
peculiar. PG 1011-040 is only slightly below the accretion rate threshold of log/L^ - -1-9 that 
we find here and, most importanly, is severely under-luminous in X-rays with respect to its multi- 
wavelength l uminosity, so that th e accretion rate derived from its 2-10 keV lumi nosity by simply 



applying the iMarconi et al.l (20041) correction, is likely severely under-estimated ([Gallagher et al 



200 ll : IVasudevan and Fabianl 120071) . NGC 7213 is the only bright Seyfert 1 galaxy known to un 
ambiguously lack reproce ssing features from Compton-thick distant material (|Bianchi et al.ll2003L 
20041 : Lobban et alJboioh . 



The inclusion of Compton-thick sources in our sample, makes the HBLR versus non-HBLR 
separation, in terms of either Eddington ratios or bolometric luminosities, less significant. Several 
non-HBLR Compton-thick sources lie now in the HBLR L Bol - A Edd plane, i.e. at accretion rate 
higher than the threshold we find for Compton-thin sources. This mixing could be, at least partially, 
due to the difficulties in assessing the intrinsic 2-10 keV luminosity of Compton-thick sources, 
because of the unknown geometry of the absorbers and reflectors @ 

However, such uncertainties should in principle affect equally over- and under-estimates of 
the luminosities. Instead, while no Compton-thick HBLR is still found below the Eddington rate 
threshold, a significant fraction (64%) of Compton-thick non-HBLR has accretion rates higher than 
this limit. Moreover, uncertainties in the exact geometry of absorbers and reflectors, could only be 
responsible for modest (factor of 2-3) over-estimates of the intrinsic bolometric luminosity and so 
of the Eddington accretion rate. This could perhaps explain the presence of few non-HBLR closely 
above the Eddington ratio threshold, but can hardly explain order-of-magnitude over-estimations 
of the bolometric luminosity, as inferred from the 2-10 keV spectra. 



If the scenario proposed by iNicastrol (|2000f) is correct, non-HBLR Compton-thick sources 
with Eddington accretion rate estimates much above the threshold, are likely to be peculiar ob- 
jects, where BLR line emission cannot be seen, not even in polarized light. These source should 
possess a BLR, but something prevents us from observing it. Diffe rent, alternative scenarios, to 
the intrinsic lack of a BLR, have already been proposed in the past ( Lumsden and Alexanderll2~o"oT 



3 An updated catalogue was used for this search. 

4 The same uncertainties are probably also responsible for the lack of HBLR and non-HBLR separation when 
indirect estimators of the nuclear bolometric luminosity (e.g. [<9lll]) are used. 
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Gu and Huan g 2002; iMartocchia and Mattl 12002 iTranl 120031 : lMoranll2007h . This could be either 



because of the lack of an unobscured population of hot electrons scattering the line emission along 
our line of sight, or because the orientation of these sources is such that our line of sight intercept 
larger portions of absorbing gas, and this covers both the nuclear so urce and at least part of the 
line -reprocessing electron region (e.g. IShu et al.ll2007l : IWu et al.ll201 ll and references therein). 



6. Conclusions 



The main results of this paper can be summarized as follows: 



• we presented evidence suggesting that accretion rate is the main ingredient which drives the 
presence of HBLRs in Seyfert 2 galaxies. By selecting a sample of 39 type 2 AGN with good 
quality spectropolarimetric and X-ray data, and for which we derived homogenous estimates 
of the mass of the central supermassive BH, we found a clear separation between Compton- 
thin HBLR and non-HBLR sources, both in luminosity (logL bo i = 43.90) and in accretion 
rate (log A edd = -1.9). A statistically similar separation is seen when bolometric luminosities 
are derived from infrared [Oiv] emission lines. 



Our results agree with the ones discussed in lNicastro et al.1 (120031) but have higher stat istical 



significance (due to the larger sample we use here, compared to lNicastro et al.1 (|2003)). and 
probably more robust due to the more accurate estimates of the BH masses that we derive 
here; 

the inclusion of Compton-thick sources, washes out the separation between HBLR and non- 
HBLR, but still no HBLR source falls below the Eddington ratio threshold. We propose that 
the presence of a significant fraction (64%) of Compton-thick non-HBLRs at accretion rates 
higher than the threshold found for Compton-thin sources, is not due to the lack of BLRs, 
but to heavy line-of-sight absorption preventing us from observing not only the direct line 
emission, but also their polarized light scattered off an, at least partly, obscured population 
of hot electrons. 
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Fig. 1 . — 2-10 keV luminosities derived from our x-ray analysis are plot against [Om] luminosities. 
A correction factor of 70 is used in the bottom panel. 
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Fig. 2. — The bolometric luminosity inferred from the 2-10 keV luminosity against the Eddington 
ratio for all the Compton-thin sources of our sample is plot in the top panel while in the bottom plot 
Compton-thick sources are introduced, using a correction factor of 70. The dotted lines represent 
the maximal separation values between the cumulative distribution of the two samples with respect 
to L Boi and A edd . 
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Fig. 3. — Stellar velocity dispersions obtained from the Call and Mg b absorption lines against 
those derived through the FHWM of the [Om] emission line (on the Y axis). The two estimates 
are very different and not correlated. Only the BH masses estimated with the first method will be 
used in this work. 
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Fig. 4. — In the top panel bolometric luminosities derived from the [Om] and Eddington ratios for 
all the 39 Compton-thin and Compton-thick sources of our sample are shown. In the bottom panel 
bolometric luminosities derived from [OlV] and Eddington ratios for 32 sources are shown. We 
used a dispersion value of cr = 0.3 dex for uncer tainties in the estim ates of the [Om] bolometric 
luminosities, based on a recent study presented in lRisaliti et al.l (|201 1|) where the flux of the [Om] 
line is used as a reliable indicator of the bolometric emi ssion of quasars. For [OlV] bolometric 
luminosities we used a dispersion value of cr = 0.35 dex (|Goulding et al.ll2010|) . The dotted lines 
represent the maximal separation values between the cumulative distribution of the two samples 
with respect to L Bo i and A ec id- 



